SPECTRAL PROPERTIES OF OPERATORS 
USING TRIDIAGONALISATION 

MOURAD E.H. ISMAIL AND ERIK KOELINK 

Abstract. A general scheme for tridiagonalising differential, difference or g-diffcrcnce op- 
erators using orthogonal polynomials is described. From the tridiagonal form the spectral 
decomposition can be described in terms of the orthogonality measure of generally different 
orthogonal polynomials. Three examples are worked out: (1) related to Jacobi and Wilson 
polynomials for a second order differential operator, (2) related to little q- Jacobi polynomials 
and Askey- Wilson polynomials for a bounded second order g-difference operator, (3) related 
to little q- Jacobi polynomials for an unbounded second order q-difference operator. In case 
(1) a link with the Jacobi function transform is established, for which we give a g-analogue 
using example (2). 



I. Introduction 

Bochner's 1929 theorem classifies polynomials which are eigenfunctions to a second order 
differential operator, see e.g. [12] where also a historical discussion and extensions to other op- 
erators can be found. Instead of looking for differential operators diagonalised by polynomials, 
we can also look for differential operators which are tridiagonalised by polynomials. Naturally, 
there is no need to restrict to differential operators, and we can also consider difference and 
g-difference operators. Neither is there a need to restrict to a differential or difference oper- 
ator of second-order. In [13] we have discussed a general framework as well as the examples 
of the Schrodinger operator with Morse potential and the Lame equation. In this paper we 
give a general approach for finding self-adjoint differential (or difference) operators which are 
tridiagonal in a suitable basis of orthogonal polynomials. Since the spectral decomposition of 
tridiagonal operators can be expressed completely in terms of, in general different, orthogonal 
polynomials, we obtain the spectral decomposition of the differential or difference operator 
involved. 

In Section [2] we sketch the general framework, and we specialise to the case where the 
basis in which the operator is tridiagonal is given by orthogonal polynomials. Since we start 
with orthogonal polynomials which are eigenfunctions to an explicit differential or difference 
operator, it is natural to start from the families in the (g-)Askey scheme. In Section [3] we 
work out the case of the Jacobi polynomials, which yields an operator which is, after a change 
of variable, for a special choice the differential operator for the Jacobi functions. Its spectrum 
can be given explicitly since the tridiagonal operator corresponds to Wilson polynomials. In 
Section H] we discuss a g-analogue of this situation involving the little q- Jacobi polynomials. 
However, due to lack of symmetry in the g-case we have to consider two cases. It turns 
out that one case can be solved within the g-Askey scheme using Askey- Wilson polynomials, 
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whereas the other case leads to orthogonal polynomials not in the g-Askey scheme. In Section 
E] we link to a result [T7J (3.3)] by Koornwinder which states that Jacobi polynomials can be 
mapped to Wilson polynomials by the Jacobi function transform, and we give a g-analogue of 
this result in Section H~T1 Koornwinder's result is predated by a result by Broad [5], which the 
limit case of the Whittaker transform mapping Laguerre polynomials to Meixner-Pollaczek 
polynomials, and is hugely extended by Groenevelt [10] to the top of the Askey-scheme by 
showing that the Wilson transform maps Wilson polynomials to Wilson polynomials. 

Starting within the Askey-scheme and requiring that the tridiagonalisation can be solved 
within the Askey-scheme leaves only the cases of the Jacobi and Laguerre polynomials. As 
stated, the first case is discussed in Section El whereas the second is a similar calculation which 
can also be obtained using limit transitions from the results in Section El The other cases in 
the Askey-scheme to which this procedure can be applied are the Meixner, continuous Hahn, 
continuous dual Hahn and Wilson polynomials. The orthogonal polynomials associated to the 
tridiagonalisation are no longer in the Askey-scheme, but these are orthogonal polynomials 
associated to birth and death processes which is a general phenomenon. A further study of 
these examples is desirable, and may lead to a better understanding of this procedure. A 
similar remark applies to the orthogonal polynomials in the g-Askey scheme, see also Section 

E 

The tridiagonalisation method is also known as the J-matrix method. Tridiagonalisation for 
different operators in all kinds of bases is a much employed method in many physical problems, 
see e.g. [2J, [TTJ, [2TJ, and references in [T3]. We hope that the more general approach given 
here widens the class of operators to which the tridiagonalisation method applies. In Section 
El we discuss briefly the possible extension of this method to an even larger class of operators. 

2. General framework 

In this section we introduce the general set-up for tridiagonalising certain differential or 
difference operators. In Section 12.11 we give the general framework which is specialised to the 
case of orthogonal polynomials in Section 12.21 

2.1. General framework. Let \x and v be positive Borel measures on the real line R so that 
v is absolutely continuous with respect to \i. Let 5 = ^ be the Radon-Nikodym derivative, 
so 5 > 0. Let P be a function space that is dense in both L 2 (/j,) and L 2 {y) with respect to 
the respective topologies. 

We assume that we have an operator L: T> — >■ T> operator such that (L,T>) is symmetric 
in L 2 (fi). Moreover, we assume the existence of an orthonormal basis ($ n )neN °f L 2 (fi) of 
eigenfunctions of L, so $„ G V for all n G N and L$„ = A n $„, A n G E. Then (L, V) is 
essentially self-adjoint. 

Next we assume that there is an orthonormal basis ((/> n )^L of L 2 {y) such that (fi n G V for 
all n G N and that 

n = A n $ n + B n $„_!, A n , B n eR (2.1) 

(with the convention B = 0). Finally we assume that multiplication by the reciprocal of the 
Radon-Nikodym derivative preserves the space P, so M: T> — > V, Mf = 5~ x f. 
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The object of study is the operator T = ML: T> — > T> as an operator on L 2 [y). Note 
(T0 n ,0 m ) L 2 H = (A n ML^ n + B n ML$„_i , A m $ m + B m $ m ^) L 2 {v) 

= (A n A„$ n + B n A n _i$ n _i, A m $ m + P m $ m _i) L 2 (M) (2.2 
= A n A n 5 n+ i5 n+ i )jn + (A n A n -\-B 2 A n _i)5 n>m + A n _]A„_i-B„5 n)m+ i 

In particular, (T,T>) is a Jacobi operator on the Hilbert space L 2 {y)\ 

T(j) n = a n (f) n+1 + b n (j) n + a n _i0 n „i 
a n = A„A n S n+1 , b n = A n A 2 + A n _ x B l n 



2 , A D 2 (2-3) 



where the convention Bq = gives a_i = 0. Since a n , b n G R and (0 n )^l o is an orthonormal 
basis of L 2 {y) it follows that T with domain D(T) finite linear combinations of the orthonor- 
mal basis (0„)^ o is a symmetric densely defined operator to which standard techniques of 
tridiagonal (or Jacobi) operators can be applied, see e.g. [T], [I], [7], [13], [H], [19] and refer- 
ences given there. We assume D(T) = T> so that we can obtain the spectral decomposition of 
suitable self-adjoint extensions of (T,D). We assume that (T,T>) is essentially self-adjoint. In 
general the spectral decomposition of the closure of (T, V) depends on the possible vanishing 
of the coefficients a n and on the growth behaviour of the sequences (a n )^ =0 , (& n )^L , see e.g. 
[131 Thm. 2.7]. 



Remark 2.1. Note that we can trivially change L by adding a suitable real constant, say 
V = L + 7, and we change A n to A 7 = A„ + 7. The operator in general changes by adding 
a multiplication operator T 7 = T + jM. By comparing the coefficients a 7 = a n + jA n B n+ i, 
b^ = b n + ■y(A 2 l + B 2 ) we find the three-term recurrence 

M0 n = A n P n+1 n+1 + (A* +^)0 n + A n _ 1 5 n n _ 1 . 

The additional freedom is useful, see Sections [3] and SJ 

2.2. Specialising to orthogonal polynomials. The examples we consider are based on 
orthogonal polynomials. Let fi and v be orthogonality measures for orthogonal polynomials; 

P n (x)P m (x) dfi(x) = H n 5 nimj / p n (x)p m (x) dis(x) = h n 5 n , m . (2.4) 

Jm. 

We assume that both /1 and v correspond to a determinate moment problem, so that the 
space V of polynomials is dense in L 2 (n) and L 2 {y). We also assume that f R f(x)dfi(x) = 
j R f(x)r(x) dv(x), where r is a polynomial of degree 1, so that the Radon-Nikodym derivative 
^ = 5 = 1/r. Then we obtain, using lc(p) for the leading coefficient of a polynomial p, 

* = ^p. + mp-%^p~-i (2.5) 

lc(P n ) f/ n _i lc{p n ) 

and by taking <j) n , respectively $ n , the corresponding orthonormal polynomials to p n , respec- 
tively P n , we see that (12. ip specialises to 



A (h R (Ti A 1C R i ( \ / K lc(Pn-l) / 9 fi N 

v4 n $ n + B n Q n _i, A n = , i — , 5 n = lc(r)W— — — — . (2.oJ 

lc(P n ) V «n V iZfi-l lc(Pn) 
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We assume the existence of a symmetric self-adjoint operator L with domain T> = V on L 2 (p) 
with LP n = A n P n , and so L$ n = A n $ n , for eigenvalues A n e R. Then (L,V) is an essentially 
self-adjoint operator, since P is dense in L 2 (p), as assumed in Section [2.11 By convention 
A_! = 0. So this means that (P n ) n satisfies a bispectrality property, and we can take the 
family (P n ) n from the Askey scheme or its g-analogue. 

Considering the operator T = rL with domain T> = V on L 2 (v) shows that the coefficients 
in (12.31) are given by 



A , / x HPn) IK+I , « R n ( lc(p n ) V , . , , x 2 _^n / lc^ -i) X ' 

By switching from L and T to L 7 and T 7 , see Remark l2.lt for a suitable constant 7 we can 
assume that A 7 7^ for all n G N. 

and a 7 7^ for n > 0, and a 1 _ l = 0. So we need to solve for the orthonormal polynomials 
r n (A) satisfying 

Ar n (A) = a 7 r n (A) + 6 7 r n (A) + a^r^A), (2.8) 

We assume that the corresponding Jacobi operator is essentially self-adjoint (or equivalently 
that the orthogonal polynomials {r n } correspond to a determinate moment problem). Thus 
we have a positive Borel measure p so that 



r n (X)r m (X) dp(X) = cW (2.9) 

In particular, this implies that we assume T with domain V essentially self-adjoint. 

Remark 2.2. Polynomials associated to a birth and death process with birth-rate A n and 
death-rate p n are generated by 

-xF n (x) = \ n F n+1 (x) - (A„ + p n )F n (x) + /x n F n _i(x), F (z) = 1, F x (x) = (A + /x - ar)/A . 
with A n > for all n and fi n > for n > 1, fi > 0, see [121 Ch. 5]. Apart from sign issues the 
recursion (12.71) corresponds to a birth and death process with birth-rate A n = |A 7 ||^- ^j^fryj 

and death-rate p n = |A^_ 1 |(lc(r)) 2 jj fcn - (^p~J^) ^ or n > 0, po = 0. Assuming that the signs 
of a 7 and A 7 are independent of n e N, this allows us to draw conclusions on the support of 
the spectrum of T 7 , i.e. the orthogonality measure of the polynomials r n , see |T2l Ch. 5, 7]. In 
this case, let e = sgn(a 7 ), r\ = sgn(A 7 ). In particular, the spectrum of T is contained in [0, 00) 
if (s,rj) = (—1, 1) or (1, 1) and the spectrum of T is contained in (— 00, 0]) if (s,rj) = (1, —1) 
or (—1, —1). In the examples A 7 is bounded from above or from below, so that it is possible 
to find a 7 G R so that all eigenvalues A 7 have the same sign. In the examples we see that 
the condition on a 7 is also valid. 

Defining U: L 2 [y) — > L 2 (p), U : n 1— > r n , gives a unitary operator satisfying UT = U*M, 
where M is the multiplication (by A) operator on L 2 (p). It follows that the spectrum of T 
equals cr(T) = supp(p), and the spectrum is simple. 
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It follows that 

OO OO / \ 

F(s;A) = J2 r nWM*) = $^r n (A)^p (2.10) 

n=0 n=0 ^ n 

is a formal eigenvector for the eigenvalue A for T. The series in ( I2.10p converges in L 2 {v) for 
A a discrete mass point of p and diverges for A in the continuous part of p. In that case the 
series, which is a non-symmetric Poisson kernel, 

OO , s 

F t (x;\) = J2t n r n (\)^=^l (2.11) 

n=0 y/fln+1 



converges for \t\ < 1. 



3. Jacobi polynomials 

We work out the programme of Section 12.21 for the case of the Jacobi polynomials and the 
related hypergeometric differential operator. 

For the Jacobi polynomials P^^\x), we follow the standard notation [3J, [T2], [T5]. We 
take the measures \i and v to be the orthogonality measures for the Jacobi polynomials for 
parameters [a + 1,(3), and (a, (3) respectively. We assume a,(3 > —1. So we set P n (x) = 
Pi a+1 ^\x), p n (x) = P„\x). This gives 

2n + a + p + 1 T(n + a + p + 1) n\ 

w n 7^ (n + a + /3 + l)„ 

lc(p n ) = Z n (a) = — — : , lc(P n ) = l n (a + 1 . 

2 n nl 

Moreover, r(x) = 1 — x. Note that we could have also shifted in (3, but due to the symmetry 
Pn (x) = (—l) n Pn (—x) of the Jacobi polynomials in a and (3 it suffices to consider the 
shift in a only. 

The Jacobi polynomials are eigenfunctions of 

Lf(x) = (1 - x 2 ) f"{x) + ((3-a-l-(a + f3 + 3)x)f'(x), 
LP (.<x+m = _ n ( n + a + [3 + 2 ) p( a+1 '® 

so that A n = -n(n + a + (3 + 2). We set 7 = -(a + 5 + l)(/3 -5 + 1), so that = 
-{n + a + 5 + l)(n + /3 - 5 + 1). Hence, we study on L 2 ([-l, 1], (1 - x) a (l + x) p dx) the 
operator (1 — x)(L + 7) which is 

d 2 d 
T 7 = (l-a;)(l-a; 2 )— + (l-x)(/3-a-l-(a+/3+3)a;)— -(l-x)(a+5+l)(/3-5+l). (3.1) 
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Applying the results of Section I2~21 gives the explicit expression for the recursion coefficients 



a 



7 



6 7 



2(n + a + 8 + l)(n + P - 5 + 1) (n + 1) (n + a + 1) (n + /3 + 1) (n + a + (3 + 1) 
2n + a + /3 + 2 y (2ri + a + /3 + 1) (2n + a + /3 + 3) 

2(n + a + 5 + l)(n + /3 - 5 + l)(n + a + l)(n + a + (3 + 1) 
(2n + a + i9 + 1) (2n + a + /3 + 2) 
2n(n + + a + 5 + l)(n + /3 - 5) 



(2n + a + /3) (2n + a + (3 + 1) 

Then the recursion relation (12.71) for |T 7 is solved by the orthonormal version of the Wilson 
polynomials 

W n {fi 2 ; -(1 + a), -(1 + a ) + «J, 1(1 - a ) + /9 - 5, -(1 + a)), 

see [15], where the relation between the eigenvalue A of T and fi 2 is given by A = —2 (^-) 2 — 
2/r\ 

Theorem 3.1. Let a > —I, (3 > —I, and assume 7 = -(a + 5 + 1)((3 -5 + 1) el. The 
unbounded operator (T 7 ,P) defined by (13.1 p on L 2 ([— 1, 1], (1 — x) a (l + x)? dx) with domain 
the polynomials V is essentially self- adjoint. The spectrum of the closure T 7 is simple and 
given by 

(-00, - l(a + l) 2 ) U {-ha + l) 2 + 2(1(1 + a) + 6 + k) 2 : fceN, -(1 + a) + 5 + k < 0} 

U {--(a + l) 2 + 2(1(1 - a) + /3 - 6 + I) 2 : leN, -(1 - a) + /9 - 5 + Z < 0} 

where the first set gives the absolutely continuous spectrum and the other sets correspond to 
the discrete spectrum of the closure of T 1 . The discrete spectrum is either empty or at most 
one of these sets is non-empty. 

Remark 3.2. Note that in Theorem 13.11 5 G R or = |(/3 — a). In the second case there is 
no discrete spectrum. 

The eigenvalue equation T 7 f\ = Xfx is a second-order differential operator with regular 
singularities at —1, 1, 00. In the Riemann-Papperitz notation it is 

f-1 1 00 

V < -|(l + Q!)+iA a + 5 + 1 x 
{-/3 -I(i + a )+iA 0-6 + 1 

with the reparametrisation A = —\{a + l) 2 — 2A 2 of the spectral parameter. In case 7 = 0, 
i.e. 8 = —a — 1 or 5 = 1 + (3, We set 5 = 1 + (3, and we put T = T°, a° = a n , 6° = 0. Note 
that we have degenerate case, since a = as well as 6 = and a n 7^ for n > 1. It follows 
that T0 O = o o 0o = 0' an d we have to study the recurrence relation (12. 7p on (C^o)" 1 , i- e - we 
start at n = 1 instead of at n = 0. The recurrence relation (12. 7p starting from n = 1 can also 
be solved similarly in terms of Wilson polynomials, namely by the orthonormal version of 

W n (n 2 ; 1(1 + a), 1(1 - a), -{3 + a),0 + -(3 + a)), 
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where the relation between the eigenvalue A of T and /i 2 is given by A = —2 (^j^-) 2 — 2u 2 as 
in the derivation of Theorem 13.11 

Proposition 3.3. The operator (T,V) defined by (13.1 j) with"/ = wi£/i domain the polynomi- 
als V C L 2 ([— 1, 1], (1 — x) a (l +x) /3 c?x) is essentially self-adjoint. The spectrum of the closure 
T is simple and given by 

(-oo, --(a + l) 2 ) U {--(a + l) 2 + 2(-(l - a) + /c) 2 : fceN, -(1 - a) + k < 0} U {0} 

where the first set gives the absolutely continuous spectrum and the other sets correspond to 
the discrete spectrum of the closure ofT. 



Remark 3.4. Note that putting 5 — 1 + (or 5 = —1 — a) in the result of Theorem 13.11 
precisely gives back the statement of Proposition 13.31 The spectral decomposition U is now 
defined on (C^o)" 1 by <fi n +i l— r n, where the r n are the orthonormal Wilson polynomials with 
parameters |(1 + a), |(1 — a), |(3 + a), /3 + |(3 + a)). 

The eigenvalue equation Tf\ = \f\ is a second-order differential operator with regular 
singularities at —1, 1 and oo, hence can be solved in terms of hypergeometric functions. 
Changing variables x = 1 — 2 cosh -2 1 transforms 2T in the second-order differential operator 

d 2 d 
— + ((2/3+1) coth t - (2a + 2(3 + 3) tanh t) — 
at 1 at 

on L 2 ((0, oo),2 a+/3+2 (sinh)f:) 2/3+1 (cosht)- 2a - 2/3 - 3 rft). This is the Jacobi differential operator 
as in [IB] , where (a, (3) of [16] correspond to ((3,— a — (3 — 2). Proposition 13.31 provides an 
alternative proof of the simplicity and location of the spectrum of the Jacobi differential 
operator, see [TOJ Thms. 2.3, 2.4], but it does not give the Jacobi functions as eigenf unctions. 

Having established the link of T to the Jacobi differential operator, we know that the Jacobi 
function transform gives the spectral decomposition. In particular, applying the integral 
transformation 

/■oo 

/(A) = / 0f- a - /3 - 2) /(t)(2sinht) 2 ^ +1 (2cosht)- 2a - 2 ' 3 - 3 rft 
Jo 

for / G L 2 ((0,oo), (2sinht) 2/3+1 (2cosht)- 2a - 2/3 - 3 cft) see [IS], [H] for details, to (E2D starting 
from n = 1 shows that the Jacobi function transform of t h-» 0„ + i(1 — cosh~ 2 t) satisfies the 
same recurrence relation as the orthonormal Wilson polynomials corresponding to Remark 
13.41 Since the solution to the recurrence relation is completely determined by its starting 
value we obtain 

/■oo 

/ <j)f- a -^ 2 \t)P^(l - 2 cosh" 2 1) (2sinht) 2 ^ +1 (2cosht)- 2a - 2/3 - 3 rft = 
Jo 

^^^n(^; \(l + «), \(l - a), i(3 + «),/? + i(3 + a)) 

/■oo 

x / ^'-"-^(^^(l - 2 cosh" 2 1) (2sinht) 2/3+1 (2cosht)- 2a - 2/3 - 3 ^, 
Jo 

where S n is the squared norm of the Wilson polynomials. However, this is a special case of a 
result by Koornwinder [161 §9], [T71 (3.3)] mapping Jacobi polynomials to Wilson polynomials. 
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In our case there is a drop in the degree of the polynomial, which does not occur in [TTJ (3.3)]. 
Applying [T7J (3.3)] to the integral on the left hand side gives a Wilson polynomial of degree 
n+1 with parameters (|(a + l), + — -(a + l), /3 + |(a + 3)). Since two of the parameters 
add up to zero the first term of the hypergeometric series vanishes and the polynomial can be 
written as a linear term times a Wilson polynomial of degree n. This can also be viewed as 
a special case of the so-called 'Diophantine' properties introduced by Calogero et a!., see [61 
(3.10)] and references given there. 

Remark 3.5. The Whittaker (confluent hypergeometric) differential operator can be obtained 
as a limit case of the Jacobi differential operator, see [T7]. The operator arises from this 
approach if we start with Laguerre polynomials instead of Jacobi polynomials. 

4. Little q- Jacobi polynomials 

Next we work out the programme of Section I2T21 for the case of the little q- Jacobi polynomials 
as a g-analogue of Section [3J 

We follow the notation of [To]. For < a < q~ x , b < q^ 1 we define the little q- Jacobi 
polynomials by 

q~ n , abq n+1 



p n (x;a,b;q) = 2 fi [ ' ',q,qx 
The little q- Jacobi polynomials are orthogonal with respect to a discrete measure on g N ; 

oo 

22 w k(a, b; q)p n (q k ; a, b; q)p m {q k ] a, b; q) = 5 n>m h n (a, b; q), 

k=0 

t , \ (qb]q)k, , k , ( , v (a6g 2 ;g)oo (1 - abq) (aq) n (q,qb;q) 7 
w k {a, 6; q) = -, r— (aq) , h n (a, 6; q) - 



(4.1) 



(q;q)k ' (ag;g)oo (1 - abq 2n+1 ) (qa,qab;q) n 

The leading coefficient is given by 

l C (p B (-;a,6; g )) = (-l)V^ (w - 1) (fl6gW+1;g)w 



(aq] q) n 

Defining the second-order g-difference operator 

(L^f)(x) = ^ (f(qx) - f{x)) + ^ (fiq-'x) - /(*)) , 

'tjC Ob 

B(x) = a(bqx — 1), D(x) — x — 1 

we have 

L (a 'V(-; a,b;q) = A n (a,b;q)p n (-;a,b;q), A n (a, b; q) = q~ n {l - q n ){\ - abq n+1 ). 

4.1. Case one: shift in a. In the notation of Section [2.21 we take P n (x) = p n (x;aq,b;q), 
H n = h n (aq, 6; q), p n (x) = p n (x; a, b; q), h n = h n (a, 6; q). We assume < a < q^ 1 and b < q^ 1 . 
Then \i and v are discrete measures and fi({q k }) = Wk(aq, b; q) = q k Wk(a, b; q) = q k v({q k }), so 
that r(x) = x. Also L = L^ aq ' b \ A n = A n (aq, b; q) and T 7 is given as a second order difference 
operator on L 2 (u) by 

(Ty)(x) = aqibqx - 1) (f(qx) - f(x)) + (x - 1) (fiq^x) - f(x)) + jxf(x). (4.2) 
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In particular, T 7 is a bounded operator on L 2 (v). 

We now replace 7 = (1 + abq 2 ) — {acq + bq/c), so that 

Al = q- n (l-acq n+1 )(l-bq n+1 /c), 

and we assume c choosen so that 76I and A 7 > for all n6l This is the case when e.g. 
bq < c < 1/aq, and by the conditions on a, b this is always possible. A calculation gives the 
values 



a? 



aq(l - acq n+1 )(l - bq n+1 / ' c) (1 - q n + l )(l - aq n+1 )(l - bq n+1 )(l - abq 



b~> 

t i 



(l-abq 2n + 2 ) y (1 -abq 2n + l ){l - abq 2n + 3 ) 

(1 - acq n+1 )(l - aq n+1 )(l - abq n+1 )(l - bq n+1 /c) (43) 
(1 -abq 2n + l )(l -abq 2n+2 ) 
aq(l - q n )(l - acq n )(l - bq n )(l - bq n /c) 
+ (1 - abq 2n ){\ - abq 2n + l ) ' 

The corresponding recurrence relation (12. 8p for ■^=T' y corresponds exactly to the orthonor- 
mal Askey- Wilson polynomials with parameters (y/aq, Cy/aq, -v/J, yfaq)- Explicitly, 



r n (A) = Cp n (x; y/aq,Cy/aq, -\/-, y/aq\q) 



bq (1- abq)(q,aq,bq;q) n 

C = {acq,—;q)r 



2n+l 



c y (1 — abq 

where A = 1 + aq — 2yfaqx and using the standard notation for the Askey- Wilson polynomials 
[3], [9], [12], [T3] . So the continuous spectrum of T is the interval [1+aq— 2y/aq, l+aq+2y/aq] = 
[(1 — y/aq) 2 , (1 + y/aq) 2 ]. Discrete spectrum in the orthogonality measure can only occur if 
I c a/^?I > 1 or |-v/J| > 1, see [9], [15]. Using the explicit form of the discrete mass points 
in the orthogonality measure for the Askey- Wilson polynomials then leads to the following 
theorem. 

Theorem 4.1. Consider the operatorT" 7 defined by (14. 2 p with r y = (1+abq 2 ) — (acq + bq/c) G 1R 
and with the assumption A 7 > for all n G N. Then T 7 acting on L 2 (u) with v the discrete 
orthogonality measure (14. ip for the little q-Jacobi polynomials with < qa < 1, b < q^ 1 is a 
bounded self-adjoint operator with simple spectrum at 



[(1 - y/Eq) 2 , (1 + y/Eq) 2 } U {(1 - q- k /c)(l - acq l+k ) : \q k Cy/aq\ > 1, k G N} 

U {(1 - bq l+1 /c)(l - acq- l /b) : \q l -^\ > 1, I G N} 



Remark 4.2. Theorem 14. II is a g-analogue of Theorem 13. II However in Theorem 14.11 we can 
have that the discrete spectrum can be empty, consist of one of these sets or consists of both 
sets. The degenerate case 7 = or c = 1/aq can be considered as a g-analogue of the Jacobi 
differential operator, cf. Proposition 13.31 
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We can also obtain the spectral decomposition of T acting on L 2 (v) in a direct way. Suppose 
, : q N ->• C Sc 

equivalent to 



Vx- q N — > C satisfies T 7 y A = \y x . Put y\(q k ) = Pk ^ ; then being an eigenfunction of T is 

y/w k (a,b;q) 



aq + 1 — A 



Pit 



(A) = V(i-g fe+1 )(i-^ +1 W(A) 



c V a 

This corresponds to the orthonormal Al-Salam-Chihara polynomials with (cy/aq, \ \J~^) as 
parameters, so 

P fc ( A ) = // , 1 x =Qkfa Cy/aq, -J-\q) (4.4) 
\/{q,bq;q;) k cy a 

where 2x = (aq+l — X) / y/aq. Here we use the notation for the Al-Salam-Chihara polynomials 
IS] ; [12], [15] ■ The spectrum equals the support of the orthogonality measure a of the Al- 
Salam-Chihara polynomials Qk(x\ Cy/aq, ^\f^\q), and we obtain another proof of Theorem 

So the spectral decomposition of T 7 on L 2 [y) is given by the unitary map V : L 2 (v) — > L 2 (a) 
(Vf)(x) = (/, Ux) l 2 (u) = Y,f{q k ) Qk ^ C f*>-^ v^M) 

-E^ fc )f^- Q k {x-,c^ q \ q -\ q ) 

where A = ag + 1 — 2y/aqx, x = cos 6*. The unitarity of / can also be checked directly using 
the orthogonality relations for the Al-Salam-Chihara polynomials. So VT 1 = MV, where V 
is the multiplication by aq + 1 — 2y/aqx. 

Applying V to the three-term recurrence for T 7 gives 

(aq + l- 2y/aqx)(V(j) n )(x) = a n (V '<fr n+1 ) (x) + b n (V(j) n )(x) + a n -i(V(f) n ^ 1 )(x), 

so that V(j) n is a constant multiple of the orthonormal Askey- Wilson polynomials with param- 
eters y/aq,Cy/aq, \ \f\-, s/aq- The constant is V(f) . Explicitly, 

(Vp n (-;a,b; q))(x) = ^p n (g fc ;a,6; q)) \—Qkfacy/aq,-J-\q) 

k=o W' q > k _ cVa (46) 

{aq) n,2 (acq n+l } bq n+l /c;q) 00 , , — , — b [q . — . . 

■p n (x; y/aq, Cy/aq, -J -, yfaq\q) 



(aq; q) n [y/aqe id , y/aqe ie ; q)^ ' ' ' c\ a'' 

where x = cos9. Since V is unitary h n 5 ntm = (V(j) n (-;a,b;q)),V(p m (-;a,b;q))) L 2^ estab- 
lishes the Askey- Wilson polynomials p n (x; y/aq, Cy/aq, v/J, \fwj\q) a s orthogonal polynomials 
with respect to 

da 



(y/aqe ie , y/aqe~ ie , y/aqe ie , y/aqe- ie ; q)^ 
which is well-known [3], [9], [12], [15]. 
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It is straightforward to prove ( 14. 6 p from the explicit expression for the little g- Jacobi polyno- 
mials and the generating function [9], [T2"| (15.1.10)], JT5J (3.8.13)] for the Al-Salam-Chihara 
polynomials, this leads to the slightly more general result 

h^ n{q } (X;C ' Aq) = WM^^\ a q ,tc,td 

(4.7) 

for \t\ < 1 and x = cosO. 

In order to get the full four-parameter Askey- Wilson polynomials as orthogonal poly- 
nomials from this interpretation we observe that f n : q k i— )■ a k p n (q k ; a,b; q) and g n : q k \— > 
Oi~ k p n {.q k ] a, b; q) for yfaq < a < 1/Waq are biorthogonal funtions in L 2 {y). By ( 14. 7p we find 



, , a n {aq) n / 2 {aacq n+1 ,abq n+1 /c;q) O0 b fq 

[Vfn){x) = r— f — —- . p n {x; a^/aq, c^/aq, -W-, y/aq a\q) (4i 



and a similar expression for Vg n . Now h n 5 nm = (Vf n , Vg m )i,2(a) gives the full four-parameter 
family of Askey- Wilson polynomials p n (x; a^/aq, c^/aq, ^/aq/a-\q) as orthogonal polyno- 
mials with respect to 

da 

(ay/aqe 19 , a y /aqe~' 10 , a -1 ^faqe %d , a~ 1 y / age~' 6 '; g)^ 



So we have given yet another proof of the well-known orthogonality relations of the Askey- 
Wilson polynomials, see [3], [9], [12], [15] . 

Remark 4.3. By the results of Section [3] we may consider ( I4.8P as a g-analogue of [161 §9], 
[T71 (3.3)] mapping the (bi-) orthogonal little g- Jacobi polynomials to the full four-parameter 
family of Askey- Wilson polynomials, cf . Remark 14.21 To see the formal limit transition of the 
Al-Salam-Chihara polynomials (with c = 1/aq corresponding to 7 = 0) of ( I4.4p to the Jacobi 
functions 0^f'~° _/3 2 ' we replace e 10 = q^, fi G [0,27r/lng]. Substituting as well a = q a , 
b = q 13 , we can write the Al-Salam-Chihara polynomial of ( 14. 4p as 

Qk(\(q^ + q-^); q-^ 1+a \ q^ l+ ^ l+a) \q) = 

22 ^ttt^ ;«)0l^«ti 



1=0 



1 



^ ( « - VO. K- 1 - « + . i _ l\ = tf-«-^\t\ z = cosh" 2 t 

This is only a formal limit transition, and we refer to [18] for another limit transitions of 
orthogonal polynomials from the g-Askey scheme to Jacobi functions. Koornwinder [18] gives 
a rigorous limit transition from the continuous g-ultraspherical polynomials to Jacobi functions 
with a = (3, where the degree of the polynomials is related to the argument of the Jacobi 
function as in the formal limit described in this remark. The Jacobi function transform is 
related to the spherical Fourier transform on rank one symmetric spaces, see [16] , and for the 
case a = b = 1, and c = 1/aq this corresponds to the occurrence of the Al-Salam-Chihara 
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polynomials in the spectral decomposition of the Laplace-Beltrami operator on the quantum 
disc, see Vaksman J2DJ §1.5.4]. 



4.2. Case two: shift in b. In the notation of Section I2~2l we take P n (x) = p n (x; a, bq; q), H n = 

h n (a, bq; q), p n (x) = p n (x; a, b; q), h n = h n (a, b; q). We assume < a < q^ 1 and b < q^ 1 . Then 

fx and v are discrete measures and fJ>({q k }) = Wk{ a i bq; q) = (i-6g) w h{ a i q) = r(q k )v({q k }), 
where r(x) = Also L = L^ bq \ A n = A n (a, bq; q) we take 7 = (1 + abq 2 ) - (acq + bq/c) 

as in Section 14.11 Again we assume A 7 > for all n G N. Then T 7 is given as a second order 
difference operator on L 2 (u) by 

(th f\t \ a ( b( l 2x ~ ~ bqx) 

(ry)(x) = — — — (f( qx ) - f(x)) + 

( '' 1,(1 Hr) (f(q~ l x) ~ f{x)) + ((1 + abq 2 ) - (acq + bq/c))^f^-f(x) 



x (1 — bq) v 1 (1 — bq) 

(4.9) 

Note that T is unbounded as an operator on L 2 (u). 
Then we can calculate 



a 7 



bqy/aq (1 - acq n+1 )(l - bq n+l /c) / (1 - ^^(l - ag"+ 1 )(l - bq n+1 )(l - abq^ 1 ) 



" 1-bq (l-abq 2n + 2 ) \j (1 - a6g 2 " +1 )(l - abq 2n + 3 ) 

_ n (l - acq n+l )(l - bq n+1 /c)(l - a6g n+1 )(l - bq n+1 ) 



n ~ q (1 - bq)(l - a6g 2n+1 )(l - abq 2n+1 ) 

,2 n+2 (1 ~ acg" +1 )(l - bq n+1 /c)(l - q n )(l - of) 
q (l-bq)(l-abq 2n )(l-abq 2n + 1 ) 

(4.10) 

Note that (a 7 ) n is a bounded sequence, even though T 7 is an unbounded operator. Since the 
sequence (£$)„ is a sum of a bounded sequence and unbounded sequence tending to infinity, 
it follows that [U Thm. VIII. 4] gives that T 7 is essentially self-adjoint. However, these 
polynomials do not fit in the g-Askey-scheme. This proves the first statement of the following 
proposition. The statement on the location of the spectrum follows, since we are dealing with 
birth and death processes, see Remark 12.21 since we assume A 7 > 0. 

Proposition 4.4. The unbounded operator T defined by ( 14. 9 p on L 2 (u) with v the discrete 
orthogonality measure ( 14. ip for the little q-Jacobi polynomials with < qa < 1, b < q^ 1 is an 
unbounded self-adjoint operator with simple spectrum located in [0, 00). 

As in Section 14.11 we can try to solve the eigenvalue equation Ty\ = \y\ directly. Putting 
y\(q k ) = pfc(A) / sj Wk(a, b; q) leads to the recurrence 

(1 - bq)Xp k (X) = a k p k+1 (\) + b k p k (\) + a fc _ip fe _i(A), 

a k = -v^g- fc -5(l - bq k+2 )^l -bq^ 1 )^ - q k+1 ), (4.11) 

b k = q~ k (l - bq k+1 )(a(l - cq k+1 ) + (1 - bq k+1 /c)). 
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Using jH Thm. VIII. 4] we see that the polynomials correspond to a determinate moment 
problem. It is remarkable that the support of the orthogonality measure of (14. lip and the one 
arising from Proposition 14.41 coincide. 

5. Concluding remarks 

We comment on several possible extensions of the results described, which are outside the 
scope of this paper. 

In the context of Section [2] we can also obtain tridiagonalisations of higher order differential 
or difference operators by looking at the operator p(L)Mp(L) for p a real-valued function, 
such as a polynomial of degree d and M a multiplication by a linear function. Then T is again 
a symmetric differential or (q-) difference operator on L 2 (fi) if L is self-adjoint on L 2 (n) and 
diagonal with a suitable basis of eigenfunctions. Note that the degree of T is 2d on L 2 (/i). A 
possible link to Krall-type polynomials needs to be investigated, see e.g. [8]. 

Secondly, we can extend the relation (12. ip to allow for more terms. This leads to higher-order 
recurrences for T, which may be solved by matrix-valued orthogonal polynomials. Conversely, 
for known operators T we can obtain spectral information on matrix-valued orthogonal poly- 
nomials that occur in this way. A similar remark is valid if we look for tridiagonalisation on 
Z instead of N. 
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